Astronomy & Astrophysics manuscript no. Variations' Atmos'HD189733b 


©ESO 2013 


March 8, 2013 





Letter to the Editor 



Temporal variations in the evaporating atmosphere of the 

exoplanet HD 189733b 

A. Lecavelier des Etangs 1 ' 2 , V. Bourrier 1 ' 2 , P. J. Wheatley 3 , H. Dupuy 1 ' 2 , D. Ehrenreich 4 , A. Vidal-Madjar 1 ' 2 , 
G. Hebrard 1 - 2 , G. E. Ballester 5 , J.-M. Desert 6 , R. Ferlet 1 ' 2 , and D. K. Sing 7 

1 CNRS, UMR 7095, Institut d'astrophysique de Paris, 98 bis boulevard Arago, F-75014 Paris, France 

2 UPMC Univ. Paris 6, UMR 7095, Institut d'Astrophysique de Paris, 98 bis boulevard Arago, F-75014 Paris, France 

3 Department of Physics, University of Warwick, Coventry CV4 7AL, UK 

4 UJF-Grenoble 1 / CNRS-INSU, Institut de Planetologie et d'Astrophysique de Grenoble (IPAG) UMR 5274, Grenoble, France 

5 Lunar and Planetary Laboratory, University of Arizona, 1541 E. University Blvd., Tucson, AZ 85721-0063, USA 

6 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 

7 Astrophysics Group, School of Physics, University of Exeter, Stacker Road, Exeter EX4 4QL, UK 



(N 



3 



ABSTRACT 



Oh 



Atmospheric escape has been detected from the exoplanet HD 209458b through transit observations of the hydrogen Lyman-a line. 
Here we present spectrally resolved Lyman-ff transit observations of the exoplanet HD 189733b at two different epochs. These 
HST/STIS observations show for the first time, that there are significant temporal variations in the physical conditions of an evapo- 
rating planetary atmosphere. While atmospheric hydrogen is not detected in the first epoch observations, it is observed at the second 
epoch, producing a transit absorption depth of 14.4±3.6% between velocities of -230 to -140kms~'. Contrary to HD 209458b, these 
high velocities cannot arise from radiation pressure alone and require an additional acceleration mechanism, such as interactions with 
stellar wind protons. The observed absorption can be explained by an atmospheric escape rate of neutral hydrogen atoms of about 
Q ' 10 9 g s -1 , a stellar wind with a velocity of 190kms~' and a temperature of ~10 5 K. 

An X-ray flare from the active star seen with Swift/XRT 8 hours before the second-epoch observation supports the idea that the 
observed changes within the upper atmosphere of the planet can be caused by variations in the stellar wind properties, or by variations 
in the stellar energy input to the planetary escaping gas (or a mix of the two effects). These observations provide the first indication 
CZ2 ' of interaction between the exoplanet's atmosphere and stellar variations. 
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\ 1. Introduction 2. Observations, data analysis, and results 

» : 

Observations of the transiting extrasolar planet HD 209458b in 21 ■ Observations 
the Lyman-Q' line of atomic hydrogen (Hi) have revealed that 
\Q ', this planet is losing gas (Vidal-Madjar et al. 2003). Subsequent To address these problems, we observed two transits of 
O ■ theoretical studies indicate that atmospheric escape (so-called HD 189733b on 6 April 2010 and 7 September 2011 with 
'evaporation') arises from the intense stellar X-ray and extreme the Space Telescope Imaging Spectrograph (STIS) onboard the 
1~l '. ultraviolet energy input into the upper atmosphere (Lammer et Hubble Space Telescope (HST). The data consist of time-tagged 
al. 2003; Lecavelier des Etangs et al. 2004; Yelle 2004), leading observations obtained with the G140M grating, producing time- 
to moderate escape rates for massive hot-Jupiters, or to forma- resolved spectra from 1 195 to 1248 A with a spectral resolution 
tion of planetary remnants when strong evaporation implies a of about 20kms~ 1 at 1215.6A (Lyman-ff) with exposure times 
dramatic change in the planet mass (Lecavelier des Etangs et al. of 1800 to 2100 seconds. Between each consecutive HST or- 
2004, 2007; Charpinet et al. 201 1). bit, data acquisition is interrupted for about 3500 seconds by the 

Earth's occultation. The obtained spectra show stellar emission 

Despite the importance of evaporation on the fate of planets lines of H , Lyman _ a> si m (120 6.5 A), O v (1218.3 A) and the 

at short orbital distances, the physics of the exosphenc gas Nv doub i et (1242.8 A and 1238.8 A), 
and role of the star-planet system properties remain debated 

(Garcia Munoz 2007; Schneiter et al. 2007; Holmstrom et al. For both the 2010 and 201 1 observations, data were recorded 

2008; Lecavelier des Etangs et al. 2008; Murray-Clay et al. during two HST orbits before the transit, one orbit during the 

2009; Ben-Jaffel & Sona Hosseini 2010; Guo 2011). This is transit, and one orbit after the transit. We measured the tran- 

further compounded by the limited number of observations sit signature of the planetary atmosphere by comparing spectra 

(Vidal-Madjar et al. 2004; Ballester et al. 2007; Ehrenreich et taken during transit to those taken before and after the transit 

al. 2008; Fossati et al. 2010; Linsky et al. 2010), which include event. For each stellar line, we calculated transit light curves of 

non-spectrally resolved Lyman-Q' observations of the exoplanet the total emission flux and of the flux within given wavelength 

HD 189733b (Lecavelier des Etangs et al. 2010). ranges; the signature of the planetary atmosphere is detected as 

an excess absorption during the planet's transit. Because no at- 
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Fig. 1. Lyman-Q emission line of HD 189733 in April 2010 and September 201 1. Spectra obtained before (black) and during the 
transits (blue) are displayed as a function of radial velocity relative to the star. The double-peaked profile results from a single 
stellar emission line absorbed at the center by interstellar hydrogen, which produces a broad absorption feature from about -100 to 
+50 km s _1 . While no transit signatures are detected in 2010, two absorption regions are detected at more than 3-cr during the transit 
of 201 1; these regions are plotted by gray zones. They are seen at the top of the red wing around +80kms~ I and, most significantly, 
in the blue wing with a ~100kms~' wide absorption around -200 kms -1 . 



mospheric signature is detected in the Si m, O v, or N v lines, 
hereafter we consider only the Lyman-Q observations. 

The Lyman-o- line is the brightest stellar line in our 
STIS spectra from 1195 to 1248 A, and has a total flux of 
about 1.8xl0 -13 erg s _1 ctrT 2 (about 10 times brighter than for 
HD 209458). With the resulting high signal-to-noise ratio, there 
is no need to co-add several observations to detect the signature 
of the atmosphere, which was necessary for the STIS observa- 
tions of HD 209458b and for the observations of HD 189733b 
with the Advanced Camera for Surveys (ACS) of the HST 
(Vidal-Madjar et al. 2003; Lecavelier des Etangs et al. 2010). 
This allows for the first time a search for temporal variations in 
the physical conditions of the planetary upper atmosphere be- 
tween two observational epochs (here 17 months apart), as per- 
formed for the deeper atmosphere using emission spectroscopy 
(Grillmair et al. 2008). 

2.2. Detection of temporal variations in the evaporating 
atmosphere 

The Lyman-o- emission line from HD 189733 is spectrally re- 
solved. At the resolution of the G140M grating, the line is com- 
posed of two peaks separated by a deep absorption due to inter- 
stellar atomic hydrogen (Fig.Q}. In the raw data, the stellar emis- 
sion line is superimposed with the geo-coronal airglow emis- 
sion from the upper atmosphere of the Earth (Vidal-Madjar et 
al. 2003). This geo-coronal emission can be well estimated and 
removed in the final spectrum using the CALSTIS data pipeline 
(version 2.32 of 5 November 2010). Independent re-analysis of 
raw data using the same methodology as for STIS observations 
of HD 209458 (Desert et al. 2004) confirms that the airglow does 
not affect our measurements. Moreover, because we used a nar- 
row spectrograph slit of 0.1", the airglow contamination is lim- 
ited to the central part of the Lyman-Q line and does not con- 
taminate the line wings where the transit atmospheric signatures 
are detected (see below). The data of September 201 1 exhibit a 
notably low airglow emission level. 

The Lyman-Q observations of April 2010 do not show an 
atmospheric transit signature. In these data, the transit depth for 
the total flux of the whole Lyman-Q line is 2.9+1.4%, which 
agrees with the 2.4% transit depth of the planet body alone as 



seen from the visible to the near-infrared (Desert et al. 2009, 
201 1; Sing et al. 201 1). In addition, we see no excess absorption 
in any portion of the Lyman-Q spectral line profile. 

However, this situation strongly contrasts with the observa- 
tions of September 201 1, in which we see an excess absorption 
in the total flux of the Lyman-Q line that yields a transit depth of 
5.0+1.3%. This level is consistent with the results obtained with 
the non-resolved HST/ACS spectra of 2007-2008 (Lecavelier 
des Etangs et al. 2010). Importantly, the line profile shows two 
deep absorption regions at specific wavelength intervals during 
the transit: first in the blue part of the spectrum from about - 
230km s _1 to -140km s" 1 , and in the red part of the spectrum 
from 60 to HOkms" 1 . 

In the blue wing of the Lyman-Q spectrum of 2011, the 
most significant absorption is visible in the range of -230 to 
-140kms _1 , which gives an absorption depth of 14.4+3.6% (4-cr 
detection) corresponding to an excess absorption due to hydro- 
gen atoms of 12.3+3.6%. The false-positive probability to find 
such an excess anywhere in the whole searched range of -350 
to -50kms~' is only 7%. In the red wing of the Lyman-Q spec- 
trum of 2011, absorption is found in the range between 60 to 
HOkms -1 , which yields an absorption depth of 7.7 ± 2.7% (3-cr 
detection) and corresponds to an excess absorption due to hydro- 
gen atoms of 5.5+2.7%. The false-positive probability to find 
such an excess over the whole searched range [40 to 200 km s -1 ] 
is 39%. This high probability shows that the absorption seen 
in the red wing of the spectrum may not be real and is possi- 
bly caused by statistical noise in the data, although interestingly 
enough, a similar absorption was also observed in HD 209458b 
(Vidal-Madjar et al. 2003). 

None of the excess absorption features detected in the 
September 2011 data are seen in the April 2010 data. We con- 
clude that there are significant temporal variations of the phys- 
ical conditions within the extended exosphere of this extrasolar 
planet between these two epochs (Fig. |2j». 

2.3. Models 

In September 2011, the absorption depth and velocity range 
show that neutral hydrogen atoms are present up to very 
high altitudes at velocities exceeding the escape velocity of 
60kms~ 1 ; this unambiguously demonstrates that atmospheric 
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Fig. 2. Plot of the flux between -230 and -140 km s _I in the blue 
wing of the Lyman-a- line as a function of time relative to the 
center of the planetary transit. Vertical dashed lines show the 
beginning and end of ingress and egress of the transit. The red 
triangular symbols are for the 2010 observations, while the blue 
square symbols correspond to observations of 2011. Horizontal 
error bars centered on the symbols show the duration of the ex- 
posures in each HST orbit. The time-tagged data allow indepen- 
dent sub-exposures to be extracted within each HST orbit (not 
shown here), resulting in the same transit signal within error 
bars. The light curve of the planet's transit at optical wavelengths 
is displayed as a solid black line. The blue dashed line shows 
the calculated flux using the numerical simulation with an EUV 
ionizing flux 5 times the solar value, a stellar wind of protons 
with a temperature r~10 5 K, a velocity v~190kms _1 and den- 
sity «~3xl0 3 ctrT 3 together with an atmospheric escape rate of 
10 9 gs-'. 



gas must be escaping from HD 189733b. For HD 209458b, the 
Lyman-a excess absorption was observed in the spectral range 
between -130 to -50kms _1 , which is readily explained by the 
stellar radiation pressure accelerating hydrogen atoms up to 
-130km s" 1 (Lecavelier des Etangs et al. 2008). The case here of 
HD 189733b, which shows excess absorption at higher velocities 
between -230 and -140 kms -1 , is more challenging to explain. 
Assuming a distance of 19.3 pc and following the method de- 
scribed in Ehrenreich et al. (201 1), we estimated the interstellar 
medium absorption and calculated the Lyman-a emission line 
profile as seen from the planet. From the extinction-corrected 
line profile we estimated that radiation pressure can accelerate 
hydrogen atoms up to a radial velocity of -120km s _1 in the ex- 
osphere of this planet (below this radial velocity the stellar flux 
at the corresponding wavelength in the core of the emission line 
is sufficiently high for the radiation pressure to exceed the stellar 
gravity). Therefore, an additional acceleration mechanism be- 
yond radiation pressure is required to explain the high radial 
velocities of hydrogen measured during the transit. Charge ex- 
change with stellar wind protons can produce the observed high 
velocities (Holstrom et al. 2008; Ekenback et al. 2010). 

To investigate this possibility and interpret the observed Hi 
light curve, we developed a numerical Monte-Carlo simulation 
of the hydrogen atom dynamics. The details of the model will 
be given in a forthcoming paper (Bourrier et al. in preparation). 
In this N-body simulation, hydrogen atoms are released from 
HD 189733b's upper atmosphere and are rapidly accelerated by 
the radiation pressure up to 120km s and then to higher ve- 
locities by charge exchange with protons from the stellar wind. 
This dynamical model provides the time-dependent distribution 
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Fig. 3. Plot of the relative absorption observed in the blue wing 
of the Lyman-Q' stellar line (blue histogram) for the transit of 
September 7, 201 1. The dashed line shows the model with radi- 
ation pressure only; in this case, there is no absorption at radial 
velocities below — 120 kms -1 , resulting in a large x 1 of 22.8. 
If a stellar wind and charge exchange is considered, hydrogen 
atoms can be accelerated to the higher observed velocities. The 
model with an escape rate of 10 9 gs~' (solid line) gives ax 2 of 
13.0 for 17 degrees of freedom. 



of positions and velocities of the escaping hydrogen atoms in 
the cloud surrounding HD 189733b. With this information, we 
calculated the corresponding absorption over the stellar emis- 
sion line and the resulting transit light curve, which can be di- 
rectly compared with the observations. We find that the observa- 
tions are well-fitted with an escape rate of neutral hydrogen of 
about 10 9 g s _1 and a stellar wind with a temperature T~10 5 K, a 
density n~3xl0 3 cm -3 , and a velocity ~190kms~'. The best-fit 
model yields a% 2 of 13.0 for 17 degrees of freedom for the ab- 
sorption spectrum given in Fig. [3] The EUV flux controlling the 
neutral hydrogen ionizing timescale should be about 5 times the 
solar value to explain the moderate absorption observed after the 
transit of the planet (Fig. |2j. 

3. Swift X-ray simultaneous observations 

The evaporation of hot Jupiters is driven by the X-ray /EUV irra- 
diation of the planet by its parent star. To quantify the level of X- 
ray irradiation at the time of our observations in September 20 1 1 , 
we obtained contemporaneous observations with the X-ray tele- 
scope (XRT) of the Swift spacecraft (Gehrels et al. 2004). A 
source was detected toward HD 189733 with a mean count rate 
of 0.01 19+0.0007 s _1 ; consistent with previous XMM-Newton 
observations (Pillitteri et al. 2010), we did not detect the M- 
star binary companion, HD 189733B, but we did find evidence 
for a very weak hard X-ray source located about 13" south of 
HD 189733A. 

The Swift XRT spectrum of HD 189733 can be fitted with 
a multi-temperature optic ally- thin thermal plasma model (Mewe 
et al. 1985; Liedahl et al. 1995) that is typical of the coronal 
X-ray emission from active stars. Using a three-temperature fit 
(temperatures of 0.12, 0.46 and 4.5 keV) we found an average 
X-ray flux in the 0.3-3 keV band of 3.6xl0~ 13 erg s _1 cirr 2 , con- 
sistent with XMM-Newton observations of HD 1 89733 (Pillitteri 
et al. 2011). This flux corresponds to a planetary irradia- 
tion rate of 1.2xl0 24 ergs _1 , which could drive a mass loss 
rate of up to 1.0xl0 n gs _I (assuming 100% evaporation ef- 
ficiency and taking evaporation enhancement by tidal forcing 
into account; Lecavelier des Etangs et al. 2007; Erkaev et 
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Fig. 4. Swift XRT X-ray light curve of HD 189733 about the time of the September 7, 2011 transit. The data were binned into 
one point per snapshot visit, with typical exposure times of about 27 minutes. Vertical dashed lines show the beginning and end of 
ingress and egress of the transit. A bright flare occurred about 8 hours before the transit. The observed average count rate during 
the flare snapshot is a factor 3.6 higher than the mean for the whole light curve, indicating a peak X-ray flux (0.3-3 keV) that must 
be less than 1.3xl0~ 12 ergs~' cm" 2 . For comparison, the right panel shows the distribution for 63 epochs of Swift measurements 
(including the 16 obtained in September 201 1), covering a wide range of timescales. The flare occurring shortly before the transit is 
the highest X-ray flux of all 63 measurements. 



al. 2009). Assuming a realistic emission measure distribution 
(Sanz-Forcada et al. 2011), we can estimate the total luminos- 
ity across the X-ray /EUV band at the time of our observation to 
be 7.1xl0 28 erg s _1 , corresponding to an energy-limited evapo- 
ration rate of 4.4x10" gs _I . The X-ray irradiation is consistent 
with the estimated escape rate, which would thus require around 
1% efficiency in the conversion of input energy into mass loss 
(Ehrenreich & Desert 2011). But this is only a lower limit be- 
cause the estimated escape rate of neutral hydrogen atoms of 
about 10 9 gs" 1 does not include the escape of ionized hydrogen 
at the exobase of the atmosphere, and it is therefore a lower limit 
for the net escape from HD 189733b. 

The Swift XRT light curve of HD 189733 (Fig.© shows that 
the star exhibits significant X-ray variability, and most notably, 
a bright flare that occurred about 8 hours before the planetary 
transit. This flare could explain the observed variations in the 
extended cloud of high-velocity hydrogen atoms escaping the 
planet, because this could affect the properties of the stellar wind 
needed to accelerate the atoms to the observed radial velocities. 
Besides, the enhanced X-ray/EUV irradiation associated with 
this flare must lead to a significantly enhanced escape rate. With 
our best-model parameters, an enhanced escape rate leads to a 
more extended exospheric hydrogen cloud and thus to a stronger 
absorption after about one hour; the absorption remains at a high 
level during a typical ionization timescale, which is constrained 
from the post-transit Lyman-a observations to be about 5 hours. 
Therefore, an X-ray flare occurring 8 hours earlier is expected to 
lead to higher escape rate that is then detectable in Lyman-a^ 

4. Conclusion 

Whether they are related to the observed X-ray flare or 
not, the temporal variations in the evaporating atmosphere of 
HD 189733b are clearly detected in Lyman-a^ The variability of 
the neutral hydrogen cloud around HD 189733b can explain the 
high dispersion of absorption depth measurements in spectrally 
non-resolved Lyman-a' observations (Lecavelier des Etangs et al. 
2010); combining this with the present high signal-to-noise ratio 
spectrally resolved observations, we conclude that escape signa- 
tures are detected in about half of the total five transits observed 
in Lyman-a. More simultaneous X-ray and Lyman-a' observa- 
tions are needed to obtain a better picture of the complex rela- 
tionship between the stellar energetic input to the planet and the 
atmosphere's response to it, and to constrain theoretical models 



of a space weather event on hot-Jupiters (e.g. Cohen et al. 201 1). 
The HD 189733 system appears to be the target of choice, but 
future observations should also enlarge the diversity of stellar 
and planetary system properties to better distinguish the effects 
of the stellar-planet interactions from the intrinsic variability in 
the observed atmospheres. 
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